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Abstract 

We have computed dipole strength distributions for nickel and tin isotopes within 
the Relativistic Quasiparticle Time Blocking approximation (RQTBA). These cal- 
culations provide a good description of data, including the neutron-rich tin iso- 
topes ^^'''^'^^Sn. The resulting dipole strengths have been implemented in Hauser- 
Feshbach calculations of astrophysical neutron capture rates relevant for r-process 
nucleosynthesis studies. The RQTBA calculations show the presence of enhanced 
dipole strength at energies around the neutron threshold for neutron rich nuclei. 
The computed neutron capture rates are sensitive to the fine structure of the low 
lying dipole strength, which emphasizes the importance of a reliable knowledge of 
this excitation mode. 
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1 Introduction and motivation 

Neutron capture cross sections are one of the essential nuclear inputs for stud- 
ies of r-process nucleosynthesis. As this astrophysical process, which produces 
about half of the heavy elements in the universe, runs through nuclei with 
extreme neutron excesses, most of the needed cross sections are not experi- 
mentally known and have to be modelled. Such estimates are usually done on 
the basis of the statistical Hauser-Feshbach [T]|2]|3]|4] model and require, as one 
of the important ingredients, the modelling of the electromagnetic transition 
strength from nuclear states above the neutron threshold. Often dipole tran- 
sitions contribute dominantly to the cross sections. In the absence of experi- 
mental data, these dipole strength distributions for neutron-rich nuclei have 
been conventionally described by Lorentzians with an energy- or temperature- 
dependent width [5.)6i)7j and parameters derived from data for stable nuclei. 
Obviously deviations from the parametrized Lorentzian at energies around the 
neutron threshold can lead to significant modifications of the neutron capture 
cross sections. As noted first by Goriely f8j, such effects can arise from the 
so-called pygmy resonance which leads to an enhancement of the low-lying El 
strength, potentially in the astrophysically important energy range. In fact, 
an enhancement of the El strength at excitation energies below 10 MeV 
has been experimentally observed for neutron-rich nuclei relative to the one 
for stable nuclei. Prominent examples are the tin isotopes where the neutron- 
rich nuclei ^^^~^^^Sn exhibit a noticeable portion of the total El strength at 
< 10 MeV, in contrast to the stable isotopes ^^^'^^°Sn [9J. Theoretically 
these tendencies are well reproduced by calculations based on the Quasipar- 
ticle Random Phase Approximation (QRPA) on top of the Hartree-Fock plus 
Bardeen-Cooper-Schrieffer (HF-BCS) [T0| or of the Hartree-Fock-Bogoliubov 
(HFB) [11] approaches and by the relativistic QRPA (RQRPA) calculations 
[T2] and are usually explained as so-called pygmy resonances in which the ex- 
cess neutrons oscillate against the isospin-saturated core. Significant progress 
has been achieved very recently by extending the RQRPA by the quasiparticle- 
phonon coupling model using the Relativistic Quasiparticle Time Blocking 
Approximation (RQTBA). As it has been demonstrated in [T3], this extended 
approach reproduces the experimental El strength distributions in general 
very well, inluding the observed low-lying strength It has been also found 
that isotopic dependence of the excitation energy and the integrated strength 
of the pygmy dipole resonance agrees well with previous theoretical studies 
[Tips] . Ref. [16j has shown that once the theoretical strength distributions 
are corrected for the experimental detector response and in particular for the 
fact that Coulomb dissociation experiments can only access strength above 
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the neutron threshold, the data for ^^"Sn and ^^^Sn [T7] are nicely reproduced. 
This motivates us to use this model to calculate the El strength distributions 
for the chains of tin and nickel isotopes, including nuclei close and on the 
r-process path for which data do not exist, yet. The focus of our studies is 
on the effect of enhanced low-lying dipole strength on neutron capture cross 
sections, in particular for neutron-rich nuclei. Our calculations are based on 
the statistical model which we briefly describe in the next section. In section 
3 we present the RQTBA dipole strength functions and compare them with 
standard Lorentzian parametrizations and with the QRPA calculations of 
Ref. [TD]. The impact of different dipole strength distributions on neutron- 
capture rates for the nickel and tin isotopes is also discussed. 



2 Formalism 

The Hauser-Feshbach expression for the cross section of an (n, 7) reaction 
proceeding from the target nucleus i in the state fi with spin and parity 
vrf to a final state z/ with spin and parity vr^ in the residual nucleus m via 
a compound state with excitation energy E, spin J, and parity vr is given by 
(see also Pl3p8] ) 

ttH'^ 1 TI^T" 

= 2M,„E,„ (2Jf + 1)(2J„ + 1) + 

where Ei^n and „ are the center-of-mass energy and the reduced mass for the 
initial system. J„ = 1/2 is the neutron spin. For (n, 7) reactions the residual 
and the compound nucleus are the same. The transmission coefficients = 
T„(i?, J, vr; , Jf , vrf ), = T^(£', J, vr; J^, vr^) describe the transitions 
from the compound state - characterized by {E, J, it) - to the initial and final 
state, respectively. The sum of the transmission coefficients over all possible 
exit-channels is given by Ttot- It is possible to derive the laboratory cross 
section (the target is in its ground state, e.g. = 0) from Eq. ([T]) by summing 
over all possible final states u in the residual nucleus. Therefore it is convenient 
to introduce transmission coefficients which only depend on the compound 
quantum numbers: 

T,(E, J, tt) = ^ T,(E, J, tt; e;;, j;;, <). (2) 

V 

Generally the discrete energy spectrum is only partially known and conse- 
quently above a certain excitation energy a level density description is ap- 
plied. In that case the sum in Eq. ([2]) translates into a discrete sum plus an 
integration over the level density PI^PS] : 
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T,(E, J, vr) = ^ T,{E, J, tt; E;;, r^, 7:^^) (3) 

!/=0 

E 

where E^ is the excitation energy of the last experimentally known state, k, 
in the nucleus m. 

In the calculation of the transmission coefficient we closely follow the for- 
malism as described in [51[T^ . In particular we consider parity-dependent level 
densities for initial, compound and final states, here following [20] . The parity- 
dependent level densities were taken from Hilaire and Goriely as calculated 
within a combinatorial approach on top of HFB single particle energies [21j . 
We note that using a standard backshifted Fermi gas level density with 
a parity-dependence as introduced in [22] gives quite similar results. 

The 7-transmission coefficient for the decay to a state v involves a sum over 
all possible photon multipolarities. For a given multipolarity, XL, and photon 
energy E^ = E — E", the transmission coefficient is related with the 7-strength 
function, /xl, via 

Txl(^, J, vr, ^^ r, r) = T^l{E, E,) = 27rEf +Vxl(^, E,). (4) 



In this paper we mainly focus on El transitions, although the Ml contributions 
have been calculated as well using the single-particle approach [2n23j . Both 
multipolarities usually dominate for neutron capture reactions at astrophysical 
energies. 

In general, the strength function depends on the initial and final states. How- 
ever, in practical applications, one usually adopts the Brink-Axel hypothesis 
which assumes that the strength function is independent of the detailed struc- 
ture of the initial state and consequently depends only on E^ and not on 
E p^25|26j . Assuming that the photoabsorption process does not depend on 
the spin of the final state, it is possible to relate the cross section, (Txl, to the 
7-strength function: 



According to Eq. ([5]), the 7-strength function and the absorption cross section 
for El transitions are connected via 

aEi(^-y) = 2Ti\hcfE^hi{E^). (6) 
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In statistical model calculations one usually assumes that the El strength 
distribution can be approximated by a Lorentzian adjusted to reproduce the 
position and width of the giant dipole resonance (GDR) 



/ei (E^ 



1 NZ 



(7) 



with A^, Z and A = N + Z being the neutron, proton and nucleon numbers of 
the compound nucleus. Here x ~ 0.2 is a factor that accounts for the neutron- 
proton exchange contribution (see e.g., The parameters E-^i and F are 
obtained by adjustment to data for the giant dipole resonance (see e.g., [28]) or, 
if data are not available, from theoretical model predictions or systematics [3]. 
The Lorentzian approach is known to be inaccurate in the low-energy region 
as the low energy tail of the Lorentzian usually overestimates the strength 
at these energies. In practical applications one therefore modifies the width 
of the GDR, F, to account for these deficiencies. Various treatments have 
been developed over the past decades (for an overview see [6j and references 
therein). Here we follow ^ and define an energy-dependent width parameter 

as: T{E^) = T^E^/Eel 

We have calculated RQTBA El strength functions as outlined in [13] and used 
them in combination with Eqs. fHl[6l) to obtain transmission coefficients for the 
El transitions. In the RQTBA, excitations in even-even nuclei are determined 
by the nuclear response function R{uj) whose matrix elements obey the Bethe- 
Salpeter equation (BSE) of the following form [T3] : 



+ Rt}'kM Y.T.wS.MkMRitMkM^ (8) 

fcsfce r)" 

where indices ki run over single-particle quantum numbers including states in 
the Dirac sea and indices ri,r]',r]" numerate forward (-I-) and backward (— ) 
components in the quasiparticle space. Pairing correlations are treated in the 
BCS approximation. The quantity R^^\uj) describes free propagation of two 
quasiparticles in the mean field and matrix elements of the amplitude W read: 

ik4,k2k3 (0))^,,'. (9) 



In Eq. (j9]) is the static part of the amplitude W based on a one-meson 
exchange interaction with a non-linear self-coupling between the mesons and 
determined by the relativistic energy functional with the parameter set NL3 
[22]. $(a;) is the dynamical part of the interaction amplitude responsible for 
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the particle-phonon coupling. It has been calculated in the quasiparticle time 
blocking approximation, see [13] for details. 

To describe the observed spectrum of an excited nucleus in a weak external 
electromagnetic field P, one needs to calculate the strength function which is 
a convolution of the response function of Eq. ([8]) with matrix elements of the 
external field operator -P^^fc^- 



In the calculations a small finite imaginary part A of the energy variable 
is introduced for convenience in order to obtain a more smoothed envelope 
of the spectrum. This parameter has the meaning of an additional artificial 
width for each excitation. This width emulates effectively contributions from 
configurations which are not taken into account explicitly in our approach. 

The dipole photo-absorption cross section, 

^Ei{E,) = -^E,Sei{E,), (11) 



is expressed through the microscopically computed strength function of Eq. ( ITOi) 
and determines the dipole 7-strength function of Eq. ([6]) which is used for the 
statistical model calculations. 



3 Results 



We have calculated statistical model cross sections for (n, 7) reactions on tin 
and nickel isotopes using four different models to describe the dipole strength 
function. Our default model (model A) is the RQTBA which, as demonstrated 
in [13], reproduces experimental photo absorption cross sections very well, in- 
cluding the low-energy cross sections, which is of special interest for the discus- 
sion here. To account effectively for the residual effect of higher configurations 
beyond RQTBA, a 200 keV imaginary part has been included in the energy 
variable of Eq. (|T0|) . As the model is currently restricted to even-even nuclei, 
we will in the following only discuss neutron capture cross sections for odd-^ 
nickel and tin isotopes. Following Refs. [SfSITU] we describe the dipole strength 
by a Lorentzian form of Eq. ([7]) with the energy dependent width parameter 
as defined above. The centroid E^i and total strength are adjusted to the 
RQTBA strength function (model B). Although the width parameter is taken 
from Ref. [5D], it gives a good reproduction of the RQTBA strength function 



6 



in the GDR region. A comparison between models A and B allows to disen- 
tangle the influence of the low-energy strength on the neutron capture cross 
sections. Finally, we have performed calculations using the Lorentzian form of 
Eq. ([7]) but with the strength parameters determined as described in Ref. [3] 
(Model C) and calculations adopting the dipole strength functions of Ref. [10], 
which have been derived by a QRPA approach on top of a Hartree-Fock-BCS 
model (model D). 




[MeV] 

Fig. 1. 7 strength functions for selected tin isotopes as calculated by the RQTBA 
approach (solid; model A) [13] and the QRPA model of Goriely and Khan 
(dashed; model D). The dotted curves show the Lorentzian parametrization of the 
strength as proposed in [3] (model C) and the dash-dotted curves correspond to a 
Lorentzian with the width taken from [3], the centroid taken from the RQTBA 
results and the total strength adjusted to the RQTBA results, as explained in 
the text (model B). The dotted vertical lines indicate the (experimental) neutron 
threshold energies. 

Fig. [H shows the dipole 7 strength functions for the four different models and 
selected tin isotopes. We note that models A and D reproduce the known 
photoabsorption cross section data for tin isotopes quite well [T3l|T6] . While 
the parametrized strength function of Ref. [3] agrees nicely with the micro- 
scopic ^^^Sn strength functions, it predicts a centroid energy which is some- 
what higher for the neutron-rich tin isotope ^^^Sn than obtained in the other 
models. 

Comparing the various model predictions it is most striking that the micro- 
scopic model predicts quite noticeable nuclear structure effects which result 
in a rather strong fragmentation of the strength function. This is mainly due 
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to the quasiparticle-phonon coupling explicitly considered in the RQTBA ap- 
proach [13]. Comparing the microscopic strength functions (model A) with 
the Lorentzian approximation of model B, we note that for most nuclei the 
parametrization overestimates the microscopic strength function at low ener- 
gies, say E < 5 MeV. This indicates that the empirical energy-dependence 
of the width parameter [5] is not supported by the microscopic model. This 
conclusion is also supported by the QRPA results which are generally slightly 
smaller than the RQTBA results at the lowest energies. The standard parametriza- 
tion utilized in [3] agrees quite well with the model B parametrization. How- 
ever, with increasing neutron excess this model predicts the centroid at higher 
energies than the other models. Consequently, the 7 strength function for ^^°Sn 
in model C is smaller at low energies than in the other models. 

As the RQTBA model includes phonon couplings which are absent in the 
QRPA model it resolves more nuclear structure details. This can lead to fluc- 
tuations in the 7 strength function around the neutron threshold and below. 
An example is given by the nucleus ^^^Sn (see Fig.[T]) where the RQTBA model 
predicts a mild enhancement of the strength just at the neutron threshold 
which is absent in the QRPA model and, of course, in the parametrizations of 
models B and C. As the neutron separation energies move to lower energies 
with increasing neutron excess within a chain of isotopes, low-lying El strength 
as observed experimentally for the tin isotopes ^^°'^^^Sn and found within 
our microscopic RQTBA approach for neutron-rich nuclei can influence the 
7 strength function around the neutron threshold significantly. A prominent 
example is ^^^Sn where the RQTBA approach predicts a strong enhancement 
of the El strength around the neutron separation energy of E = 7.3 MeV. 
Such an increased strength is not predicted by the other models used in this 
paper (the relativistic QRPA calculations of Ref. [12] predict also additional 
strength at low energies) and, as we will see below, has a noticeable effect on 
the neutron capture cross sections. For ^^°Sn the RQTBA model also predicts 
a strong fragmentation of the 7 strength at low energies E > 5 MeV with 
enhanced strength compared to the strengths of the other models. However, 
for this nucleus the neutron threshold is already quite low {Eth = 3.5 MeV) 
and this enhancement will only have modest effect on the cross sections. 

By inserting the 7 strength functions into Eqs. (1) and (4) we have calcu- 
lated cross sections for neutron capture on the ground states of the isotopes 
105,115,131,139 gj^ arising solely from dipole transitions. The results are shown 
in Fig. [2] for neutron energies of up to 20 MeV (note that the maximum 
neutron energies reached in the r-process discussed below are of the order of 
80 — 100 keV). In order to understand the impact of various magnitudes of 
7 strengths it is necessary to realize what the relevant 7-ray energies are. It 
is well known that capture reactions proceeding via the compound reaction 
mechanism result in a 7-ray cascade with a dominant contribution coming 
from 7 rays with only a few MeV, below the neutron separation energy Sn- 
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En,,ab. [MeV] 

Fig. 2. Cross sections for (n, 7)-reactions on various tin isotopes based on El strength 
functions as derived in different models, see Fig. [1] and text. 

Recently, Ref. [31], assuming the validity of the Brink- Axel hypothesis and a 
Lorentzian parametrization of the dipole strength, has shown that this also 
holds for neutron-rich targets with dominant El strength leading to levels at 
3-4 MeV below Sn unless the density of levels accessible by El transition is 
low. In the latter case, 7 transitions with the maximally possible energy of 
Sn + En may dominate. This nicely explains the dependences we find in our 
calculations. For the cases of neutron captures on ^°^Sn and ^^^Sn the dom- 
inating 7 transitions display an energy of -E^ ~ 3 + MeV with dominant 
El transition endpoints at ~ 3 MeV below Sn- Scanning through En, we can 
nicely see in Figs. [T][2] how the larger strength functions of the two parametrized 
descriptions (models B and C) result in larger cross sections than those ob- 
tained with the microscopic approaches RQTBA and QRPA (models A and 
D), until the predictions converge close to 5'^. The RQTBA model predicts 
a strong fragmentation of the dipole strength at low energies. However, this 
is not seen in the cross sections because the relevant range of available 7 en- 
ergies averages and smoothes the impact of the fluctuations. The small kink 
appearing in the ^^^Sn cross section at En = 0.5 MeV is due to the opening of 
the {n,p) channel. 

Consistent with [21], we find a different behavior in the cross sections of ^^^Sn. 
The low level density of the doubly magic compound (and final) nucleus ^^^Sn 
causes the dominance of 7 transitions directly to the ground state, with E-y = 
Sn + En- Up to -E^ < 1 MeV the influence of the RQTBA pygmy peak explains 
why the RQTBA yields the by far largest cross sections up to that energy. 
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Above that energy, the RQTBA 7 strength drops below the strengths predicted 
by the other models and consequently the resulting cross section becomes 
the smallest. After another 2 — 3 MeV (note that Fig. 1 uses a linear and 
Fig. 2 a logarithmic energy scale), the RQTBA predicts a larger 7 strength 
again and consequently its respective cross section becomes the largest again. 
Interestingly, a second effect is overlaid on top of the general strength function 
behavior and giving rise to the decrease in cross sections found for all models 
around 2 MeV neutron energy. The level density description [21] predicts a 
relatively low density of negative parity states at excitation energies around 
9 MeV in ^^^Sn. In consequence, magnetic dipole transitions to the low-lying 
positive parity states dominate (not shown in Fig. 2). This is an example for 
the importance of the complete inclusion of parity-dependent level densities. 

Except for the parity effect, the behavior of the ^^^Sn(n, 7) cross section can 
be explained similarly to the one of ^^^Sn. The level density below Sn is very 
low for such a neutron-rich nucleus and therefore the relevant 7 energies are 
= Sn + En- The strong pygmy peak predicted by the RQTBA is also 
clearly seen in the cross section. Smaller fluctuations are dampened and the 
general strong decrease in the cross section at the high energy end is typical 
for Hauser-Feshbach cross sections, as in the other reactions considered here. 
In addition, for ^^^Sn we find that the Ml contribution to the capture cross 
section exceeds the dipole part at low neutron energies (say En < 100 keV). 
Thus, at these energies the total cross sections are less affected by the dif- 
ferences in the dipole strength functions. However, our statistical model code 
uses a rather simple description of Ml transitions based on a single particle 
treatment [2||23] . A more appropriate treatment, which includes the concen- 
tration of the Ml strength around E = 7-10 MeV due to spin excitations and 
in addition considers the pronounced orbital contribution seen in deformed 
nuclei (scissors mode |32j) is called for. 

The classical r-process operates at a typical stellar temperature of order 10^ K [3] 
To estimate which impact the enhanced low-energy dipole strength could have 
on neutron capture rates under astrophysical conditions, we have calculated 
the stellar neutron capture rate {av) , where a is the neutron capture cross sec- 
tion, V is the relative velocity of the fusing nuclei in the initial channel and the 
symbol (...) indicates proper averaging over the Maxwell-Boltzmann velocity 
distributions of the two nuclei at the stellar temperature and accounting for 
the thermal population of excited states. In Fig. [3] we compare ratios of neu- 
tron capture rates as obtained in the four different models at a temperature of 
T = 10^ K. The capture rates have been computed neglecting Ml transitions 
due to their rather simple treatment in our statistical code (see paragraph 
above). We observe that rates calculated with the parametrized El strengths 
are usually larger than our RQTBA results with the noticeable exception 
for the isotopes ^^°'^^^Sn (capturing the neutron on ^^^'^^^Sn). With the addi- 
tional observation that the other microscopically calculated strength functions 
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Fig. 3. Ratios of the stellar neutron capture rates at a temperature T = 10^ K for 
tin isotopes as calculated within models B, C, and D relative to those obtained in 
model A (RQTBA). All the calculations consider only El gamma transitions. For 
descriptions of the models, see text. 

(QRPA, model D) generally predicts smaller rates than the parametrization 
models one might conclude that the parametrization of the dipole strength 
functions with the energy dependence as defined above seems to predict too 
large strength at low gamma energies (except for the most neutron-rich tin 
isotope, ^^°Sn, considered in this work) and hence overestimates the rates, 
compared to microscopic models. Interestingly, we do not observe by com- 
parison of the four models that the low-lying strength is more important in 
neutron-rich nuclei as, for example, the parametrized strength of model B ex- 
ceeds the RQTBA results by about the same amount in the neutron deficient 
and stable tin isotopes as it does in the very neutron-rich ones (^^'^^^^''Sn). 



We have also performed calculations of dipole strength functions, cross sections 
and neutron capture rates at T = 10^ K for the even neutron-rich nickel iso- 
topes ^®~^^Ni. A comparison of the rates obtained in our four different models is 
shown in Fig. HI The findings are quite similar to those for the tin isotopes. For 
all isotopes the rates calculated with the parametrized dipole strength (models 
B and C) exceeds those obtained on the basis of the RQTBA strength func- 
tions, again pointing to a non-adequate description of the low-lying strength 
by the energy dependent width parameter as chosen in [3] . Except for the cap- 
ture on ^^Ni, the agreement between the two microscopic approaches (models 
A and D) is quite satisfactory. The Lorentzian parametrization slightly over- 
estimates the microscopic rates. We note, however, that for the heaviest nickel 
nuclei the parametrization used within Ref. |3] predicts the centroid of the 
giant dipole resonance at slightly higher energies. If corrected, the downwards 
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Fig. 4. Ratios of the stellar neutron capture rates at a temperature T = 10^ K for 
nickel isotopes as calculated within models B, C, and D relative to those obtained 
in model A (RQTBA). All the calculations consider only El gamma transitions. For 
descriptions of the models, see text. 



shift in energy would result in somewhat larger rates. 



In general, our calculations for the capture rates on tin and nickel isotopes 
show a rather close agreement between the models suggesting that the en- 
hanced low-lying dipole strength predicted by theoretical models and experi- 
mentally observed for selected nuclei has only modest influence on the neutron 
capture rates. The largest effect we find for the capture on the tin isotopes 
i29,i3igj-^ (caused by low-lying dipole strength in i30'i32gj^-j where our micro- 
scopic RQTBA model predicts noticeable enhanced strength just around the 
energies of the neutron threshold, leading to an increase in the capture rate 
by factors 2-5 compared to the other models considered here. 

The coincidence that the enhanced El strength lies around the neutron thresh- 
old in ^^^Sn and a large portion of the transitions goes to the daughter ground 
state (due to suppression of excited states at modest excitation energies in 
this double-magic nucleus) results in a large enhancement of the reaction rate 
in the RQTBA model compared to the other models. As the RQTBA strength 
function for ^'^^Sn nicely agrees with data (see Fig. 1 of Ref. \iQl), this change 
in the capture rate should be adopted in future r-process simulations. 
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4 Conclusions 



We have obtained El strengtli functions for nickel and tin isotopes based 
on the microscopic RQTBA approach which has been proven previously to 
give a good description of data, including those for the neutron-rich tin iso- 
topes ^^°'^^^Sn. These data show El strengths at energies around the neutron 
threshold which are enhanced compared to the expectations from a Lorentzian 
parametrization of the giant resonance. As it has been pointed out before that 
such enhanced dipole strengths could have noticeable effects on neutron cap- 
ture rates for r-process nuclei, we have calculated the respective capture cross 
sections for the relevant nickel and tin isotopes (®^~^''Ni and ^^^~^^^Sn) and 
have compared the RQTBA results with rates obtained either within the non- 
relativistic QRPA on the basis of Hartree-Fock plus BCS calculations or by 
empirical parametrizations. Usually we find rather good agreement between 
the two microscopic models. Noticeable exceptions are ^^^'^^^Sn, and in partic- 
ular the capture rates to ^^^Sn are noticeably larger within the RQTBA model 
than in all other models, which is caused by the enhanced dipole strength in 
this nucleus at energies above the neutron threshold. We note, however, that 
our calculation of the capture rates are based on the Brink-Axel hypothesis, 
which is common in statistical model evaluations of neutron capture rates for 
r-process applications. Within this hypothesis an enhanced dipole strength 
around the neutron threshold translates to enhanced neutron capture cross 
sections only if transitions to the ground state dominate. If the capture leads 
to excited states instead the strength function at appropriate lower energies is 
relevant. However, the transition strength for excited states could be different 
than for the ground state at the moderate energies relevant to neutron capture, 
invalidating the Brink- Axel hypothesis. Consequently, tests of this hypothesis 
for dipole transitions, similar to those performed for Gamow- Teller transitions 
in Ref. [SB] are welcome. 

Two final remarks are in order. At first, the RQTBA model predicts fragmen- 
tation and low-lying strengths also in the neutron-rich isotopes, however, not 
at the low energies close to the neutron threshold. Whether further correla- 
tions beyond the RQTBA model will further push strength to lower energies 
and hence might enhance the dipole cross sections at low neutron energies is an 
open question and must wait until appropriate nuclear models can be applied 
to such nuclei. Secondly, we have also calculated the contributions of other 
multipoles to the neutron capture rate and find that for the tin isotopes be- 
yond ^^^Sn, which have quite low neutron thresholds. Ml transitions dominate 
in our calculations making the observed differences in the dipole cross sections 
unimportant for the total cross sections due to our use of parity-dependent 
level densities \2U\. However, we note that our statistical model code describes 
Ml transitions in the so-called single-particle treatment. The calculation of the 
Ml transitions should generally be improved by using Ml strength functions 
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which resemble the experimentally observed energy distribution with strong 
spin contributions around excitation energies of order 7-10 MeV and orbital 
contributions at low energies, including the pronounced scissors mode [32j in 
deformed nuclei at excitation energies of E ^ 2-3 MeV, i.e. around neutron 
threshold energies in very neutron-rich nuclei. 
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